Superconducting Analogues of Quantum Optical Phenomena: Macroscopic Quantum 
Superpositions and Squeezing in a SQUID Ring. 



o 
o 

(N 

o 

Q 

(N 



M. J. Everitt, 1 T. D. Clark, 1 P. B. Stiffcll, 1 A. Vourdas, 2 J. F. Ralph, 3 R. J. Prance, 1 and H. Prance 1 

1 Centre for Physical Electronics and Quantum Technology, 
School of Science and Technology, University of Sussex, Brighton, Sussex, BNl 9QT, U.K^\ 
2 Department of Computing, Bradford University, Bradford, West Yorkshire, BD7 1DP, UK. 
Department of Electrical and Electronic Engineering, 
Liverpool University, Brownlow Hill, Liverpool, L69 3GJ, U.K. 
(Dated: February 1, 2008) 

In this paper we explore the quantum behaviour of a SQUID ring which has a significant Josephson 
coupling energy. We show that that the eigenfunctions of the Hamiltonian for the ring can be used 
to create macroscopic quantum superposition states of the ring. We also show that the ring potential 
may be utilised to squeeze coherent states. With the SQUID ring as a strong contender as a device for 
manipulating quantum information, such properties may be of great utility in the future. However, 
as with all candidate systems for quantum technologies, decoherence is a fundamental problem. In 
this paper we apply an open systems approach to model the effect of coupling a quantum mechanical 
SQUID ring to a thermal bath. We use this model to demonstrate the manner in which decoherence 
affects the quantum states of the ring. 

PACS numbers: 74.50+r 85.25.Dq 03.65.-w 42.50.Dv 
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Introduction 

In two recent publications 0, Q we reported on the 
theoretical description of a quantum mechanical SQUID 
ring (here, a thick superconducting ring enclosing a single 
Josephson weak link device) coupled to quantised elec- 
tromagnetic field (em) oscillator modes. In this work we 
emphasised that the SQUID ring could be used to control 
various quantum phenomena involving each of the circuit 
components of the coupled system via the static magnetic 
bias flux (<f> x ) applied to the ring. These included fre- 
quency conversion between the em modes and quantum 
entanglement extending across the system, both with 
relevance to emerging quantum techn olog ies based on 
Josephson devices HH IE IE S H, 0, EH E3, E3| 
Furthermore, work by Friedman et al on SQUID rings 
has highlighted another phenomenon of potentially great 
significance to these incipient technologies, namely the 
creation of externally controlled superpositions of macro- 
scopically distinct states in a SQUID ring, or other, 
Josephson weak link based, circuit configurations [l4j . 
As will become apparent, the creation and control of such 
states is a natural application for a SQUID ring. 

In this paper we consider the creation and control of 
macroscopic quantum superposition states in a SQUID 
ring alone, uncoupled to any em oscillator modes. First, 
we consider the spectral properties of the ring Hamil- 
tonian. Then we observe that at certain points in the 
bias flux ($ x ) applied to the ring the eigenfunctions of 
this Hamiltonian form macroscopic superposition states. 
We show that a strong enough level of dissipation may 
destroy the quantum nature of these states, whilst leav- 
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ing the flux in the SQUID ring in a statistical mixture 
of two macroscopically distinguishable states. Following 
this we demonstrate that a SQUID ring can be used to 
create (form) a controllable macroscopic superposition of 
states. In addition we show that a SQUID ring with a 
sufficiently large Josephson coupling term in its potential 
can be used to squeeze coherent states. In this it is ap- 
parent that physical phenomena associated with SQUID 
rings, and with quantum circuits built around SQUID 
rings, have analogies with effects well known in the field of 
quantum optics. Indeed the SQUID ring can be viewed as 
a non-linear medium which, for example, can be utilised 
to generate entanglements, frequency conversion, super- 
position states and squeezing. However, the SQUID ring 
has significant advantages over the generally weakly poly- 
nomial non-linear media of quantum optics which are 
usually weakly coupled to external electromagnetic (em) 
fields. Thus, it is extremely non-perturbative in nature 
(and concomitantly capable of inducing extremely non- 
linear behaviour [J, 0, uM ) with a coupling to cm modes 
that can be adjusted by means of an external bias flux. 
This would appear to make the SQUID a prime candi- 
date for future developments in what is, in effect, highly 
non-perturbative quantum optics, albeit at much lower 
frequencies. In practice these frequencies would typically 
be much less than a THz for low critical temperature 
superconductors. 



Background 
Wigner and Weyl Functions 

Although the Wigner and Weyl functions are fa- 
miliar to those working in the field of quantum op- 
tics ^E 03 E3l > their use in the quantum description of 
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FIG. 1: Illustrative example of (a) Wigner and (b) absolute value of the Weyl functions of the first four energy eigenstates of 
the simple harmonic oscillator (increasing in energy from left to right). 



Joscphson weak link circuits, and in paqrticular SQUID 
rings, appears to be rather limited. The Wigner function 
is defined to be 



W{x,p) = 1-J dc(x+^( 
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where p is the density operator describing the state of the 
system with conjugate variables position (x) and momen- 
tum (p). Physically, the Wigner function can, to some 
extent, be considered as a generalisation of the wavefunc- 
tion of the quantum system under study in which we are 
provided with information in both position and momen- 
tum space. We note that the Wigner function may, and 
often does, take on negative as well as positive values. 
An important and characteristic feature of the Wigner 
function is that the quantum correlations between the 
macroscopically distinct components of a macroscopic su- 
perposition state can be seen in an obvious and graphi- 
cal way, i.e. these correlations will appear in the Wigner 
function as interference terms between the states of the 
superposition in the x — p phase plane. 

By contrast, the Weyl function is defined as 
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It is apparent here that the Weyl function of a state is 



equal to the overlap of the displaced state with the orig- 
inal state so that X and P are considered as increments 
in position and momentum. As can be seen, the Weyl 
function is a generalised autocorrelation function; it is 
also the two dimensional Fourier transform of the Wigner 
function. Thus, just as the Wigner function highlights 
the regions in the x — p plane where the wavefunction 
amplitude is significant, the Weyl function tells us where 
there exits a significant amplitude for correlations be- 
tween intervals of Ax (= X) and Ap (= P) in this plane. 
A more detailed discussion and review of Wigner and 
Weyl functions, and the relationship between them, can 
be found in the literature [lgl . 

For those unfamiliar with these functions we provide 
a specific example in figure ^ Here, we have plotted the 
Wigner functions for the first four energy eigenstates of 
the simple harmonic oscillator together with the abso- 
lute values of their associated Weyl functions. We note 
that the Weyl function is, in general, complex valued. In 
this paper, therefore, we only ever plot its absolute value 
since, for our purposes, this provides us with sufficient 
information about the correlations of the wavefunction. 



The SQUID Ring Hamiltonian 

Over the last two decades SQUID rings, viewed as sin- 
gle, macroscopic, quantum objects, have been the subject 
of considerable attention theoretically. In early studies 
the focus was primarily on time independent properties 
and the interaction of SQUID rings with external envi- 
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ronments 0, I20I l2lj . Of late there has been much inter- 
est in time dependent behaviour, for example, in solving 
the time dependent Schrodinger equation (TDSE) for a 
SQUID ring in the presence of a microwave field |22T.l23T|. 
Recently, significant efforts have been devoted to the 
experimental measurement and control of macroscopic 
quantum superposition states in SQUID rings |14| . In 
this paper we proceed to develop a theoretical descrip- 
tion of macroscopic quantum superpositions in SQUID 
rings, borrowing on techniques that are commonly used 
in quantum optics. We extend the usefulness of this de- 
scription by considering quantum mechanical squeezed 
states in SQUID rings. For both superposition of states 
and squeezing in these rings we also discuss the effect of 
dissipation (decoherence). 

In the widely used lumped component model of a 
SQUID ring yjj |2(| the Hamiltonian takes the form. 



H 



q 2 (s-^r 

2C " 



2A 



hv cos 



2tt$ 



(1) 



where $ and Q are, respectively, the magnetic flux 
threading the ring and the electric displacement flux 
across the weak link (with [$, Q) = ih), hv/2 is the ma- 
trix clement for Joscphson tunnelling through the weak 
link (with critical current I c = lev), $0 = h/2e is the 
superconducting flux quantum and A and C are, respec- 
tively, the ring inductance and the capacitance of the 
weak link in the ring. 

Introducing a unitary translation operator T = 
exp (—i<$> x Q/h), we can then write down the ring Hamil- 
tonian as 



h' = r f m 
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where it is clear that as hv — > the system behaviour re- 
duces to that of a simple harmonic oscillator. Given the 
relation between our system and the simple harmonic os- 
cillator we now define creation and annihilation operators 
in the usual way, i.e. as 



$ _ 



These raising and lowering operators, as used in quantum 
optics, then allow us to write the ring Hamiltonian in a 
more convenient form. We also choose to express it in 
dimensionless units, normalised to Hu>, where uj/2t: = 
I/271VAC is the SQUID ring oscillator frequency. This 
takes the form 









v 






■\) 


cos 








LU 





2Cu 



2TTip. z 



(3) 

where ip x = <S> x /<&q is the normalised static bias flux ap- 
plied to the SQUID ring. We note that the cosine term in 
the Hamiltonian introduces non-linearities to all orders. 
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FIG. 2: Potential energy in units of Hu) \ lu = 1/vACJ 

dimensionless flux (see equation [1J for a SQUID ring with 
parameter values C = 5 x 10~ 15 F, A = 3 x 10" 10 H and 
hv = 0.047$o/A for $ x = 0,0.49 and 0.5$ - Also shown are 
the probability density functions of the rings wavefunctions 
displaced by their energy eigenvalues. 



We have seen that this property of the SQUID ring intro- 
duces highly non-perturbative effects 0, |3] when coupled 
to other circuit systems. In this paper we show that it 
also gives rise to quantum superpositions of macroscopi- 
cally distinct states and squeezing within the ring itself. 

From © it is apparent that as long as the ratio v : w, 

and the product Clo \ = ^/cjh^j , remain the same the 

physics of this system is unchanged. We therefore choose 
values of C and hv (or cquivalently I c = 2ev) that can 
be attained using currently available micro-fabrication 
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techniques, that are physically sensible and that will 
lead to SQUID ring systems exhibiting quantum be- 
haviour at experimentally accessible temperatures. With 
these factors in mind we choose the circuit parameters 
C = 5 x 1CT 15 F, A = 3 x 1CT 10 H and a sufficiently 
large value of hv (= 0.047$o/A; I c = 2/iA) to generate 
clear wells in the ring potential. Thus, for a thin film 
Joscphson tunnel junction weak link with a lnm oxide 
insulator thickness (dielectric constant wlO) a capaci- 
tance of 5 x 10 -15 F yields junction dimensions « 0.25^m 
square, readily achieved using micro-fabrication. Again, 
with these dimensions the supercurrent density in the 
junction is around 4kAcm~ 2 which is perfectly reason- 
able. Furthermore, with C = 5xl(T 15 F, A = 3xlCT 10 H, 
u)/2tt = 130GHz, well below the frequency corresponding 
to the superconducting energy gap in niobium (w ITHz), 
a metal often used in weak link device fabrication. Given 
these chosen parameter values, and assuming, as our ex- 
ample, SQUID circuits based on niobium, these corre- 
spond to v/uj = 7.9 and Cu> = 4.1 x 10 -3 , values which, 
unless otherwise stated, we now keep to throughout this 
paper. 

Adopting these values of A, C and hv we 
show in figure the SQUID ring potential 
U($,$ x ) = {<$> -<S> x f /2k - 7iz/cos(27r$/$o) - see 
JTJ - computed for three different values of & x 
(= 0$ (top), 0.49$ (middle) and 0.5$ (bottom)). 
We also show in this figure the probability densities 
of the wavefunctions of the ring as a function of <f> x . 
These probability densities are displaced by their energy 
eigenvalues, found by solving the time independent 
Schrodinger equation (TISE) using the Hamiltonian 
|QJ. As the bottom of each of the local wells in the 
ring potential in figure [21 is approximately quadratic, 
we would expect the solutions deep within these wells 
to look like those for the simple harmonic oscillator. 
In addition, we find that, on average, these states are 
slightly squeezed in terms of the magnetic flux variable 
$. For example, in figure El the lowest state in each of 



the wells for $ x = 0.49$ has (a y/Cuj/m 



0.43 



and A 



« 0.58 (compared to 0.5 in 
the unsqueezed state), where $ and Q are normalised 

1 

to (H/Cuj) 2 and (HCuj) 2 , respectively. To aid in the 
description of the system we adopt the term locally 
s-harmonic to describe this behaviour, where the prefix 
s denotes the squeezed nature of these states. The 
equivalence of the low lying set of energy eigenvalues 
deep in each (local) well is a sound approximation except 
where the energy levels of two or more wells align. In 
such cases, of course, symmetric and antisymmetric 
superpositions of the cigenfunctions for the isolated wells 
develop. 

In figure|21we show a set of computed SQUID ring en- 
ergy eigenvalues of the Hamiltonian operator ©, start- 
ing with the ground state, as a function of external flux 
<f> x and spanning one flux quantum. For clarity the en- 
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FIG. 3: Computed SQUID ring energy eigenvalues of the 
Hamiltonian for the SQUID ring of figure |5] 



ergy levels are shown as alternating black and grey lines. 
As can be seen, for the values of A, C and hv we have 
adopted in this work there exist many crossing points 
of the original eigenvalues as $> x is changed over a $o 
period. At the scale provided in figure |3 these crossing 
points appear to be degenerate in energy (e.g. in the low- 
est two eigenvalues at $> x = 0.5$o)- Of course, this is not 
the case, as would be evident if these crossing points were 
computed to sufficient accuracy. However, for the poten- 
tial wells shown in figure |21 with the very weak coupling 
between levels in different wells, the energy splittings at 
these points may be extremely small indeed. 



Results 



Quantum superpositions of macroscopically distinct 
states 



Ignoring, for the moment, the special cases where the 
locally s-harmonic oscillator states are degenerate (zero 
coupling between wells) we instead consider making a 
equal superposition of two energy eigenstates. We start 
by assuming that these two states are locally s-harmonic 
in different wells and, for simplicity, take these states 
to correspond to the lowest energy levels in each well. 
Then, as one would expect, we create a superposition of 
macroscopically distinct states in flux. As this superpo- 
sition state is no longer an eigenstate of the Hamiltonian 
for the system, it must evolve with time. This evolu- 
tion introduces a phase difference between the two sta- 
tionary (locally s-harmonic) states in the superposition 
which manifests itself in a time dependent evolution of 
the interference term in the Wigner function. 

In computing Wigner functions for the SQUID ring it is 
convenient to introduce equivalent dimensionless position 
(x) and momentum (p) operators in place of $ and Q. 



These are defined as 

An example of a macroscopic superposition state in a 
SQUID ring, as illustrated through its computed Wigner 
function, is shown in figure 01 both in perspective and in 
projection on the x — p plane. Here, we have taken our 
standard values of A, C and hv (above) and have selected 
the state (\A) + \B)) at the flux bias $ x = 0.49$ of 

figure Gib) ■ In this example we can distinguish in the 
Wigner function two macroscopically distinct flux states 
of the SQUID ring in the x—p plane (i.e. the gaussian-likc 
components), separated from one another by an oscilla- 
tory region. The latter, oscillatory, region arises because 
of the quantum coherence between the two separate com- 
ponents in the superposition and demonstrates that we 
are indeed dealing with a true quantum superposition 
of macroscopically distinct states in flux. We note that 
the Wigner function of figure 0] has been calculated at a 
fixed time t = 0. However, its general form does not vary 
with time. Nevertheless, there is dynamical evolution of 
the interference term in the superposition but not in the 
observable flux states. 



Tunable superposition of states in a SQUID ring 

We now consider the potential with a static bias flux 
$ x = 0.5$o, as shown in the bottom plot of figure [3 By 
examining the first two energy eigenstates we find that 
we have wavefunctions which are a symmetric, |s), and 
antisymmetric, \a), superposition of the lowest energy, 
locally s-harmonic, oscillator vacuum states of the two 
middle wells. These two states, which are squeezed, are 
extremely close in energy due to the height of the barrier 
between the wells in the potential. Even so, the inter- 
ference terms in the Wigner function between the states 
for each of these macroscopically distinct states are tt/2 
out of phase. However, from the viewpoint of theory we 
might also consider an equal superposition of these states 
of the form 

-L(| S )+exp(i0)|a» (5) 

where 8 is the phase. Let E s and E a be the eigenenergies 
corresponding to the eigenstates \s) and |a) so that AE = 
E s — E a the energy splitting associated with the lifting of 
the degeneracy concomitant with tunnelling between the 
adjacent wells. Then 8 = AEt. It is seen that this phase 
8 changes as a function of time and any particular value 
of exp(i6*) is realised cyclically with period 2ir/ AE. If we 
wished to fix our state in a given superposition we can 
do this simply by changing our external magnetic bias 
flux by a fraction. If this fractional change in bias flux is 
sufficient this will remove the coupling between the wells 
and set 8. The effect of changing 8 can be dramatic, as 
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FIG. 4: Wigner function showing, both in perspective and in 
projection on to the x—p plane, the nature of this quantum su- 
perposition of macroscopically distinct states (\A) + \B)) 
for the SQUID ring of figure|21 The dimensionless quantities 
x and p are defined in equation (0J 




FIG. 5: Wigner function of a superposition of the lowest two 
(symmetric and anti-symmetric) energy eigenstates as a func- 
tion of phase - see equation © and related text and figure 0] 
for scalings. 

illustrated in figure [5] where the Wigner function in the 
x — p plane is shown for selected values of this parameter. 
As is apparent, when the wavefunction for the ring is 
strongly (but not completely) localised in two or more 
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regions in the ring potential it is the quantum interaction 
between these regions that is responsible for the creation 
of superposition of macroscopically distinct states. 

We note that given the very slightly different energies 
between the symmetric and antisymmetric superposition 
states, due to the height of the barrier between the two 
wells, the superposition (J5J will oscillate slowly back and 
forth between these wells. Using these energies we cal- 
culate this period to be 100ns. This is much longer than 
the time constant corresponding to the AC frequency of 
the SQUID ring in our example (i.e. for C = 5 x 10~ 15 F, 
A = 3x 10 -10 H this is 7.6ps) but well within the decoher- 
cncc times of modern SQUID ring circuits 0, 0,113,123 • 

In the above discussion we have considered the devel- 
opment of macroscopic superposition states in a SQUID 
ring based on a choice of ring parameters which can 
be realised by fabrication and which are physically rea- 
sonable. However, this choice does not connect di- 
rectly with published experimental data. We will there- 
fore deviate briefly from our standard parameter val- 
ues and consider an explicit example of superpositions 
of SQUID ring states as reported recently in the liter- 
ature by Friedman et al [Tj|. In this paper the authors 
considered a "Quantum Superposition of Distinct Macro- 
scopic States" ', and presented experimental evidence in- 
dicating that a SQUID ring could be placed into a su- 
perposition of two magnetic flux states. We now demon- 
strate by computation that this superposition may in- 
deed form a true macroscopic superposition state. The 
experimental system used a SQUID ring with the cir- 
cuit parameters C = 1.03 x 10~ 13 F, A = 2.38 x 10" 10 H 
and I c = 2.02 x 10~ 6 A. To obtain a superposition state 
as demonstrated in this paper we used the external bias 
flux quoted by the authors, i.e. & x = 0.514466$ . The 
potential energy of a SQUID ring with these parameters 
is shown in figure together with the probability den- 
sities for the ring wavefunctions. These are displaced, 
as before, by their energy eigenvalues. In this figure we 
have marked with an arrow the states in the left hand 
and right hand wells from which we will form our super- 
positions. These arc the states that were utilised in the 
Friedman et al experiment. 

These cigenfunctions are similar to those of figure [S] 
in that they form a symmetric, |s), and antisymmetric, 
|a), superposition of the of the locally s-harmonic states 
of the separate wells, albeit with higher and different or- 
dinal numbers. The quantum state of the SQUID ring, 
as reported by Friedman et al. 01 w hl thus be a super- 
position of these two eigenstates. For our purposes it is 
therefore sufficient to look (again) at superpositions of 
the form 

-L(| s ) + exp(;0)|a)). (6) 

In figure we show the Wigner functions for three dif- 
ferent superpositions of these eigenstates, again both in 
perspective and projection on to the x — p plane. We 
notice that compared with the Wigner functions of fig- 
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FIG. 6: Potential energy in units of hu) versus dimensionless 
flux (see equation of a SQUID ring with parameter values 
C = 1.03 x KT 13 F, A = 2.38 x 1(T 10 H, I c = 2.02 x 10~ e A and 
<fr x = 0.514466$o (after Friedmann et al [3). Also shown are 
the probability density functions of the rings wavefunctions 
displaced by their energy eigenvalues. The arrow indicates 
the states used to calculate the superposition state. 

ure[5]both the states of the macroscopic quantum super- 
position corresponding to the left and right hand wells 
in figured an d the region of interference between them, 
display more complex patterns. Given the choice of more 
highly excited states in the Friedmann et al experiment 
this is to be expected. Nevertheless, figure demon- 
strates that the Wigner function (and, of course, the as- 
sociated Weyl function) can expose sophisticated quan- 
tum coherent behaviour in SQUID rings. 

Dissipation 

In considering the effect of dissipation on the calcu- 
lations presented in this paper we have chosen to use a 
standard approach, and one familiar in quantum optics. 
This is to done by coupling the SQUID ring to a decoher- 
ing monochromatic thermal bath. The master equation 
for the evolution of the density operator of the system 
then takes the form [27| , 

%=-\\^iP\ + ^(M + l)(2opo t -o t op-/w t a) 
+ 2h M ( 2at(0a - aflt/9 _ /9aat ) • ( 7 ) 

where Al is related to the temperature T and the 
frequency tub of each decohering bath via Mi = 
(cxp (hiOb/kBT) — 1) and ji is the coupling (damping 
rate) between each of the components to its respective 
thermal bath. For the following examples that we will 
now calculate we set T = IK with ujb = 

To illustrate the effect that dissipation has upon a 
macroscopic superposition state in a SQUID ring, we 
now solve the master equation (JJJ for the system evo- 
lution, where the ring is taken to be in its ground state 
at $2, = 0.5$o- For this computation we return to our 
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FIG. 7: Wigner function, from left to right, for the three distinct phases 0,7r/2 and tt (cf. figure and after Friedmann 
et al showing, both in perspective and in projection on to the x — p plane, the superposition of the (symmetric and 

anti-symmetric) energy eigenstates as indicated in figure 



initial circuit values of C = 5 x 1CT 15 F, A = 3 x lfr 10 H 
and hv (= 0.047 < f>g/A) . Here we have chosen to use a 
decoherence rate of O.Olw. In figure [S] we have plotted 
the Wigner and Wcyl function for the ring at particular 
times. As the system evolves we notice that the Wigner 
and Weyl functions display very clearly the disappear- 
ance of the quantum coherence between the two states 
of this superposition. We also note that after sufficient 
time has elapsed the Wigner function still displays two 
distinct flux probabilities. However, the SQUID ring is 
longer in a pure state and this represents a classical coin 
toss probability and not a quantum one. This is also 
clearly reflected in the disappearance of the symmetri- 
cally positioned correlation peaks in the Weyl function 
at these later times. 



Squeezed states of a SQUID ring 

Just as in quantum optics, where coherent light can be 
squeezed in the number-phase plane through its interac- 
tion with a non-linear optical medium, so we shall now 
demonstrate that a SQUID ring, with the cosine Hamilto- 
nian l[T|). can be used to squeeze an initial coherent state. 
This starting condition in a SQUID ring can be achieved 
by changing the weak link structure in the ring. As is well 
know, the combined Josephson critical current in a par- 
allel, two weak link, loop , connected by superconducting 
wires (a DC SQUID j2£|), can be varied by adjusting the 
magnetic flux threading the loop. Such a structure can 
serve as the adjustable weak link within a larger diam- 
eter SQUID ring. This has been already been discussed 
in several earlier papers, including that of Friedman et al 
[Tij . The arrangement is depicted in figure HJ Provided 
the control magnetic flux ($ l2 ) threading the minor loop 
in figure is large enough, and the weak links in the 
loop can be fabricated (in principle) to possess identi- 



cal critical currents, then the net supercurrent carried 
through the loop can be made vanishingly small. In this 
situation the Josephson coupling energy in Q reduces 
essentially to zero leaving just the parabolic background 
term (i.e. as for a simple harmonic oscillator). Let us 
now consider a possible mechanism by which we could 
create a coherent state in a SQUID ring. If we allow the 
SQUID ring to decohere to its ground state, then, as we 
are in the low Josephson coupling limit, this state will 
very much resemble that of the vacuum state of the sim- 
ple harmonic oscillator. Now, a coherent state is just the 
vacuum state displaced from the origin and then allowed 
to evolve freely. We can realise this for a SQUID ring 
by very rapidly changing the external bias flux so that 
the wavef unction does not change significantly in shape, 
but the relative expectation value of flux (position) no 
longer corresponds to the minimum of the potential well. 
We can, within reason, choose the ramp rate, start and 
end bias flux, to create a state which is a very good ap- 
proximation to a coherent state at a required energy. We 
take this as the initial prepared condition of the major 
SQUID ring (with a main control flux $> Xl ) for which we 
can choose a coherent state \A = i\/T). Wc follow this 
initial set up by a very rapid reduction in the $ I2 to 
yield the desired (and finite) Josephson coupling energy 
in the principal SQUID ring corresponding to the ring 
potential considered throughout this paper. Clearly, in 
an experimental situation this reduction must take place 
on a time scale short compared with the decoherence time 
of the SQUID rin g, co upled to its environment. From the 
literature [T3. 124 , l25l |2(| it appears that the coherence 
times in SQUID rings can be sufficiently long to make 
this readily achievable. With this proviso the system is 
prepared in the required coherent state and then allowed 
to evolve over time using the master equation JJ}. 

Here, as before, wc take the SQUID ring circuit param- 
eters to be C = 5 x 1CT 15 F and A = 3 x 1(T 10 H. How- 
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FIG. 8: Wigner (left) and Weyl (right) functions for the 
ground state of the SQUID ring of evolving in the pres- 
ence of dissipation (T = IK and 7 = O.OIcj). The effect of 
dissipation is to remove the quantum correlations in this su- 
perposition of macroscopically distinct states over the time 
shown. 



ever, in order to improve the effects of squeezing by the 
ring we have increased hv to 0.24<E>g/A, yielding a super- 
conducting critical current density of around 20kAcm~ 2 . 
From our given initial condition, and by computing the 




FIG. 9: Schematic of a SQUID ring with a (flux) controllable 
Josephson tunnelling energy. Here the area enclosed by the 
outer loop dSi is much greater than that enclosed by the inner 
loop dS2- 



uncertainties in flux, we see that the SQUID ring can 
squeeze the magnetic flux within the ring. Wc choose 
not to show the results for the charge as they do not 
provide any additional information relevant to this dis- 
cussion. In figure 1101 we show the uncertainty in flux 
for a selection of decoherence rates. Here, the horizontal 
line denotes that squeezing happens for (Ax) 2 less than 
1/2. We can clearly see that, on average, the ring has 
effectively squeezed this state. We also see that for a ju- 
dicious choice of decoherence rate and interval of time, 
dissipation can assist in the squeezing process. We note 
that in the 7 = 0.1 case it can be seen that the SQUID 
ring decoheres to its (squeezed) vacuum state over a rela- 
tively short period. Thus, it appears that these potential 
wells may be used with facility to generate squeezing of 
the the magnetic flux variable of the ring. It is therefore 
reasonable to assume that adiabatic changes of this po- 
tential would, unless matching conditions exist between 
local s-harmonic energy levels of adjacent wells, allow 
some adjustment of the expectation value of the flux in 
the ring. 



Conclusions 

In this paper we have explored two applications of 
SQUID rings, with strong analogies to the field of quan- 
tum optics, which may prove to be of great utility given 
the current interest in quantum circuit technologies. In 
the first we describe the manner in which we can cre- 
ate macroscopic superposition states in SQUID rings and 
the way these may be manipulated through an exter- 
nal control flux. In the second we provide a dynamical 
mechanism for inducing squeezing of the magnetic flux 
variable in a SQUID ring starting from an almost har- 
monic oscillator state. These take advantage of the highly 
non-perturbative quantum nature of SQUID rings arising 
from the cosine coupling energy term in the ring Hamil- 
tonian. As we emphasise, both are of great interest in the 
light of current research in the fields of quantum comput- 
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ing and quantum information processing 0, IToL ITTl IT3| . 

Moreover, the results presented in this paper emphasise 
that the SQUID ring can act as a versatile quantum 
device and, as we have shown previously, can be used 
to create correlations across extended , m ulti-component, 
quantum circuit systems 0, d [2iJ lM HH H3| • These cor- 
relations can, for example, be made manifest as quantum 
entanglements and quantum frequency up/down conver- 
sion [J , with each controlled by the external bias flux 
applied to the SQUID ring. In this sense these phe- 
nomena, including macroscopic superposition states and 
squeezing, highlight the role of the SQUID ring as the es- 
sential machinery for a range of applications in supercon- 
ducting quantum circuit technologies. This further en- 
hances our understanding of these non-pertubative quan- 
tum objects and their possible usefulness in these incipi- 
ent technologies. 
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FIG. 10: Uncertainty in dimensionless flux (see equation 
versus dimensionless time (ut), for a SQUID ring with pa- 
rameter values C = 5 x 10" 15 F, A = 3 x 10~ 10 H and 
Hv — 0.24<l>o/A> computed for a selection of decoherence rates 
at T = IK. Here the horizontal line denotes that squeezing 
occurs for(Aa;) 2 less than 1/2. The grey background oscil- 
lation for the 7 = 0.001, 0.01 and 0.1 examples is the 7 = 
pattern repeated for comparison. 
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